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1. Introduction 

The lidar equation is a special form of the radiative transport 
equation in single scattering approximation and describes the return 
signal of a lidar. Based on this lidar equation satellite backscatter 
lidar returns P have been simulated 1 using realistic optical parameters of 
the atmosphere (ft = total backscatter coefficient, 6"= total extinction 
coefficient). The lidar equation reads 

V('r) ~ Q ' 4z. ' (~ Z! J 

where r is the range from the lidar to the scattering volume, and Cq is 

the instrumental constant including the pulse energy. Because of the 

viewing mode vertically down to the earth surface the range r is at that 

range where the return signal is detectable and the maximum range r is 
the satellite altitude. m 

The simulated lidar signal shows the bounds of atmospheric layers 
with different optical density with a height resolution of about 100 m. It 
is even possible to detect the upper bound of the planetary boundary layer 
(PBL) through a thin cirrus cloud. 'Subvisible' clouds not detectable with 
passive remote sensing techniques will be recognized with satellite lidar. 

The inversion of satellite lidar data to the profile of the extinc- 
tion coefficient and so to the optical depth of each atmospheric layer and 
the transmission of the atmosphere contributes e.g. to studies on the 
radiation flux and the radiation budget, from which heating/cooling rates 
can be derived. 


2. Sensitivity study on Klett's algorithm 

Klett's 2 algorithm for retrieving the total extinction coefficient 
profile has been developed for application to ground-based lidar returns 
by a backward integration from the far end r m to the near end r Q , the 
range where the incident and backscattered pulse overlap totally. Our 
study results in an assessment of the applicability of Klett's algorithm 
to satellite backscatter lidar returns. We used the simulated data 1 of a 1 
J Alexandrite laser operated at about 0.7 pm and at a satellite flight 
level of 840 km. Klett assumes a power law for the (i- fT-relation 
(fi = c-fr k ) , the exponent k has been chosen to be 0.7 according to Fenn 2 . 
Klett's solution with backward integration mode reads 
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while the forward integration mode reads 

^ ( (S('r)-S M)/k ) 
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where 6(r m ) is the total extinction coefficient at the far end bound (just 
above the earth surface) and 6"(ro)is the near end bound value (top of the 
atmosphere) for a satellite lidar system. 

Variations of these boundary values result in a parallel shift of 
the retrieved extinction profile ( 1 in-log-scale) to greater values for 
greater 61^) or £Jtr 0 ). Higher values of the exponent k makes the re- 
trieved profile steeper. The stability of both solutions is moderate, 
because in the forward integration mode the nominator increases and the 
denominator decreases, while the reverse is true in the backward inte- 
gration mode. 


3. Results 

In Fig. 1 the vertical profiles of the total extinction coefficient 
,3" scale ) versus the height (lin-scale) are shown resulting from Eqs.(2) 
and (3) with correct boundary values of_6(r m ) and £Ftr 0 ), respectively. 

The solid lines are the known input ^^-profiles of two atmospheric 
models. In an atmosphere with low turbidity (Fig. la) the forward inte- 
gration mode leads to a better fiRyt-profile (dashed) than the backward 
integration mode (dashed-dotted 6^ U fP r o f ile) • On the other hand in an 
atmosphere with high turbidity (Fig. 1b) the backward integration leads to 
better results. But these two examples give exceptionally good inversion 
results . 

Both modes have been applied to a large number of simulated satel- 
lite signals of multiply layered atmospheres, which differ by the aerosol 
type or the optical depth of each layer. Typical results of the inversions 
are given in Fig. 2. Fig. 2a) shows the eTj n - and two fo u t-P rofil, 5 S for 
a clear atmosphere with a cirrus layer (optical depth = 0.3) between 7.5 
and 9.5 km height, Fig. 2b) shows the set of extinction profiles for an 
atmosphere with low turbidity in the stratosphere and troposphere but very 
high turbidity in the PBL (desert dust), and Fig. 2c) shows the set of ' 

profiles for a turbid stratosphere (volcanic aerosol) and a clear tropo- 
sphere. All examples demonstrate that both inversion algorithms are 
unsatisfying for satellite lidar signals, because the inverted profiles 
are incorrect to about one order of magnitude. A reason for these unsatis- 
fying results is that the values of the pair (c,k) of the power law are 
different for different aerosol types and that c is assumed to be ccnstant 
in the derivation of Eqs . ( 2) and (3) and only one value of k is used in 
the inversion algorithm. 

Recently Klett 4 proposed an improved algorithm distinguishing 
between air molecules (subscript R) and aerosol particles (subscript M). 

The lidar equation for these two scatterers reads 
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A linear law for the B-6^relation, the aerosol lidar ratio SM(r) = 
®M(r)/flM(r) is assumed in the solution to the aerosol extinction coeffi- 
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cient. The 1 i da r ratio may change from layer to layer according to the 
aerosol type. The aerosol extinction coefficient profiles 


6 - - A(*) + 


WO 


j -h 

with A ( r ) = S M (r)-ftp(r) and X ( r) = S^( r) -P( r) -r 2 and 


X(-r)e*p(J T Jam - g~ e m) °l-t) 

+ 2{/V-e* P (fAU- 


)((«)■ e *P \ Jr (AW)- 6 r (")) dr 


6 -„W) = -AW) + 


AM7 -Z jxn ^p(-QAW)-5^))ol-r)^ 


n ■& ’u 

are the solutions of the backward and forward integration modes of the 
1 idar equation (4) . 

Above the PBL the values of A(r) and the second term (called B(r)) 


are at least one order of magnitude greater than the value of BmIt). 

Slight errors in both terms lead to great errors in 6^. If B(r) is less 
than A(r), negative values of d m( r) are gained, which is more likely with 
the backward integration mode. Fig. 3 gives the (TM-profiles according to 
Eq.(5) (dashed-dotted) and according to Eq . ( 6 ) (dashed) and the correct gT. 
input profile (solid), the simulated P(r) profile is without noise, the M 
boundary values are the correct ones, the lidar ratio S|v|(r) used is by a 
factor 0.5 too little. The figure shows a relative deviation of the re- 
trieved profile from the input profile to about 0.43 at the surface for 
the forward and to about 6.7 at the top of the atmosphere for the backward 
integration mode. The problem how to get the correct values of the lidar 
ratio has not been solved up to now; perhaps it can be estimated from the 
lidar signal itself, where the increased information of a multiple wave- 
length lidar would certainly help. 

It is not always advantageous, as Fernald- 3 states, to apply the last 
calculated G^(r) as a new boundary value 5~M( r m) or 6TM( r o) because, if 
this ^M-value is wrong, this error will propagate through the whole 
inversion. 

In fact, the problem of inverting the lidar returns may be even more 
complex because there is some evidence® that even in optically rather 
thin atmospheres the second order scattering cannot be ignored, so that an 
extended version of the lidar equation must be used. Then there is cer- 
tainly no chance anymore to get an analytical solution of the inversion 


problem. 

So we try to solve the lidar equation to the aerosol extinction 
coefficient with an iterative method, which seems to be more stable with 


respect to inaccuracies in the input data. 
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